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Abstract. Endeavours of the unification of the four fun- 
damental interactions have resulted in a development of 
theories having cosmological solutions in which low-energy 
limits of fundamental physical constants vary with time. 
The validity of such theoretical models should be checked 
by comparison of the theoretical predictions with ob- 
servational and experimental bounds on possible time- 
dependences of the fundamental constants. 

Based on high-resolution measurements of quasar spec- 
tra, we obtain the following direct limits on the average 
rate of the cosmological time variation of the fine-structure 
constant a: 

| a/a\ < 1.9 x 10~ 14 yr" 1 

is the most likely limit, and 

| a/a\ < 3.1 x 10~ 14 yr" 1 

is the most conservative limit. 

Analogous estimates published previously, as well as 
other contemporary tests for possible variations of a 
(those based on the "Oklo phenomenon" , on the primor- 
dial nucleosynthesis models, and others) are discussed and 
compared with the present upper limit. We argue that the 
present result is the most conservative one. 

Key words: quasars: absorption lines - cosmology: the- 
ory 



1. Introduction 

The fine-structure constant a = e 2 /Tic is the key param- 
eter of Quantum Electrodynamics (QED). Initially it was 
introduced by Sommerfeld in 1916 for describing the fine 
structure of atomic levels and corresponding resonance 
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lines. Later it became clear that a is important for de- 
scription of the gross structure of atomic and molecular 
spectra as well as fine one. Moreover, now we know that 
any electromagnetic phenomena may be described in terms 
of the powers of a. In reality, a is not a true constant. It 
is established in the quantum field theory and confirmed 
by high-energy experiments that the coupling constants 
depend on distance (or momentum, or energy) because of 
vacuum polarization (see, e.g., Okun 1996). Here we con- 
sider the low-energy limit of a, namely, its variation in the 
course of the cosmological evolution of the Universe. The 
possibility of fundamental constants to vary arose in the 
discussion of Milne (1937) and Dirac (1937). (For a sketch 
of the history of the problem of variability of the funda- 
mental constants, see Varshalovich & Potekhin 1995). The 
most important Dirac's statement was that the constancy 
of the fundamental physical constants should be checked in 
an experiment. 

The value of the fine-structure constant is known with 
rather high accuracy ~ 4 x 10~ 9 . The CODATA recom- 
mended value based on the 1997 adjustment of the fun- 
damental constants of physics and chemistry is equal to 
a = 1/137.03599993(52). However, even this high pre- 
cision does not exclude the possibility that the a value 
was different in early cosmological epochs. Moreover, some 
contemporary theories allow a to be different in different 
points of the space-time. Endeavours of unification of all 
fundamental interactions lead to a development of mul- 
tidimensional theories like Kaluza-Klein and superstring 
ones which predict not only energy dependence of the con- 
stants but also dependence of their low-energy limits on 
cosmological time. 

Superstring theories. The superstring theory is a real 
candidate for the theory which is able to unify gravity with 
all other interactions. At present this theory is the only 
one which treats gravity in a way consistent with quan- 
tum mechanics. In the low energy limit {E <C -Epianck = 
y / hc 5 /G ~ 1.2 x 10 19 GeV, where G is the gravitational 
constant, h is the Planck constant, and c is the speed 
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of light), the superstring theory reduces to the classical 
General Relativity (GR) with, however, an important dif- 
ference: All versions of the theory predict the existence of 
the dilaton which is a scalar partner of the tensorial gravi- 
ton. This immediately leads to an important conclusion 
that all the coupling constants and masses of elementary 
particles, being dependent on the dilaton scalar field </>, 
should be space and time dependent. Thus, the existence 
of a weakly coupled massless dilaton entails small, but 
non-zero, observable consequences such as Jordan-Brans- 
Dicke-type deviations from GR and cosmological varia- 
tions of the fine structure constant and other gauge cou- 
pling constants. The relative rate of the a variation is 
given by 

-^kH^-C^m) 2 (1) 

a 

where k is the main parameter that determines the effi- 
ciency of the cosmological relaxation of the dilaton field 
towards its extreme value 4> mi and H is the Hubble con- 
stant (Damour & Polyakov 1994). In principle, the varia- 
tion of a defined by Eq. (1) depends on cosmological evo- 
lution of the dilaton field and may be non-monotonous as 
well as different in different space-time regions. 

Kaluza-Klein theories. In superstring models the tran- 
sition from (4 + £>)-dimensional string objects to the four- 
dimensional observed reality proceed through a compacti- 
fication of the extra dimensions. Generalised Kaluza-Klein 
theories offer another possibility of unification of gravity 
and the other fundamental gauge interactions via their 
geometrization in a (4 + D)-dimensional curved space. 
In such theories the truly fundamental constants are de- 
fined in 4 + D dimensions and any cosmological evolu- 
tion of the extra D dimensions would result in a variation 
of the fundamental constants measured in the observed 
four-dimensional world. For example, in the Kaluza-Klein 
theories the fine-structure constant a evolves as 

a cx R~ 2 (2) 

where R is a geometric scaling factor which characterises 
the curvature of the additional Z?-dimensional subspace 
(Chodos & Detweiler 1980; Freund 1982; Marciano 1984). 

There are many different versions of the theories de- 
scribed above and they predict various time-dependences 
of the fundamental constants. Thus, bounds on the varia- 
tion rates of the fundamental constants may serve as an 
important tool for checking the validity of different theo- 
retical models of the Grand Unification and cosmological 
models related to them. 

In the next section we briefly consider the basic meth- 
ods allowing one to obtain restrictions on possible varia- 
tions of fundamental constants. In Sect. 3, an astrophys- 
ical method is considered in more detail. In Sect. 4, we 
specify the properties of observational data required to 
reach a high accuracy of evaluation of a/ a. In Sect. 5, we 



describe an observational technique employed to obtain 
data of the required quality. An upper bound on \a/a\ is 
presented and its theoretical consequences are discussed 
in Sect. 6. In Sect. 7, conclusions are given. 

2. Characteristic features of methods for checking 
possible variations of the fundamental constants 

Theoretical and experimental techniques used to investi- 
gate time variation of the fundamental constants may be 
divided into extragalactic and local methods. The latter 
ones include astronomical methods related to the Galaxy 
and the Solar system, geophysical methods, and labora- 
tory measurements. 

An interest in the problem of changing constants has 
increased after an announcement of a relative frequency 
drift observed by several independent research groups us- 
ing long term comparisons of the Cs frequency standard 
with the frequencies of H- and Hg + -masers. Such drift, in 
principle, could arise from a time variation of a because H- 
mascr, Cs, and Hg + clocks have a different dependence on 
a via relativistic contributions of order (Za) 2 . A detailed 
description of the laboratory method based on such com- 
parisons may be found in Prestage et al. (1995), Breakiron 
(1993). 

Stringent limits to a variation have been presented by 
Shlyakhter (1976) and Damour & Dyson (1996), who have 
analysed the isotope ratio 149 Sm/ 147 Sm produced by the 
natural uranium fission reactor that operated about 2x 10 9 
yrs ago in the ore body of the Oklo site in Gabon, West 
Africa. This ratio turned out to be considerably lower than 
that in the natural samarium, which is believed to have 
occurred due to the neutron capture by 1 Sm during the 
uranium fission. Shlyakhter (1976) has concluded that the 
neutron capture cross section in 149 Sm has not changed 
significantly in the 2 x 10 9 yrs. However, the rate of the 
neutron capture reaction is sensitive to the position of 
the resonance level E r , which depends on the strong and 
electromagnetic interactions. At variable a and invariable 
constant of the strong interaction (that is just a model 
assumption), the shift of the resonance level would be de- 
termined by changing the difference of the Coulomb en- 
ergies between the ground state 149 Sm and the excited 
state of 150 Sm* corresponding to the resonance level E r . 
Unfortunately, there are no experimental data for the ex- 
cited level of 150 Sm* in question. Damour & Dyson (1996) 
have assumed that the Coulomb energy difference between 
the nuclear states in question is not less than one be- 
tween the ground states of 149 Sm and 150 Sm. The latter 
energy difference has been estimated from isotope shifts 
and equals ~ 1 MeV. However, it looks unnatural that a 
weakly bound neutron, captured by a 149 Sm nucleus to 
form the highly excited state 150 Sm*, can so strongly af- 
fect the Coulomb energy. Moreover, the data of isomer 
shift measurements of different nuclei indicate that the 
mean-square radii (and therefore the Coulomb energy) of 
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the charge-distribution for excited states of heavy nuclei 
may be not only larger but also considerably smaller than 
the corresponding radius for the ground states (Kalvius & 
Shcnoy 1974). This indicates the possibility of violation of 
the basic assumption involved by Damour & Dyson (1996) 
in the analysis of the Oklo phenomenon, and therefore this 
method may possess a lower actual sensibility. 

Another possibility of studying the effect of changing 
fundamental constants is to use the standard model of the 
primordial nucleosynthesis. The amount of 4 He produced 
in the Big Bang is mainly determined by the neutron- 
to-proton number ratio at the freezing-out of n<->p reac- 
tions. The freezing-out temperature T/ is determined by 
the competition between the expansion rate of the Uni- 
verse and the /3-decay rate. A comparison of the observed 
primordial helium mass fraction, Y p = 0.24 ± 0.01, with a 
theoretical value allows to obtain restrictions on the differ- 
ence between the neutron and proton masses at the epoch 
of the nucleosynthesis and, through it, to estimate relative 
variation of the curvature radius R of the extra dimensions 
in multidimensional Kaluza Klein-like theories as well as 
the a value (Kolb et al. 1986; Barrow 1987). However one 
can notice that different coupling constants might change 
simultaneously For example, increasing the constant of 
the weak interactions Gf would cause a weak freeze-out 
at a lower temperature, hence a decrease in the primor- 
dial 4 He abundance. This process would compete with the 
one described above, therefore, it reduces sensibility of the 
estimates. Finally, the restrictions would be different for 
different cosmological models since the expansion rate of 
the Universe depends on the cosmological constant A. 

A number of other methods are based on the stellar 
and planetary models. The radii of the planets and stars 
and the reaction rates in them are influenced by values 
of the fundamental constants, that offers a possibility to 
check the variability of the constants by studying, for ex- 
ample, lunar and Earth's secular accelerations, which has 
been done using satellite data, tidal records, and ancient 
eclipses. Another possibility is offered by analysing the 
data on binary pulsars and the luminosity of faint stars. 
A variety of such methods has been critically reviewed by 
Sisterna & Vucetich (1990). All of them have relatively 
low sensibility. 

An analysis of natural long-lived a- and [3- decay- 
ers, such as 187 Rc or 149 Sm in geological minerals and 
meteorites, is much more sensitive (e.g., Dyson 1972). 
A combination of these methods, reconsidered by Sis- 
terna & Vucetich (1990), have yielded an estimate a/a = 
(-1.3 ±6.5) x 10~ 16 yr" 1 . 

The following weak points are inherent to the methods 
described above: 

(i) The derived restrictions strongly depend on model con- 
jectures. 

(ii) The "local" methods give estimates for only a nar- 
row space-time region around the Solar system. For 
example, the time of the operation of the Oklo reac- 
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Fig. 1. Fine Doublet Splitting from Quasar Spectra 



tor corresponds to the cosmological redshift z k> 0.1, 
while a method based on observations of quasar spec- 
tra (described below) makes it possible to bring out 
evolutionary effects for large span of redshifts toz~ 5. 

This and other problems, as we believe, can be solved 
by an astrophysical method which will be described in the 
next section. 



3. The astrophysical method 

An astrophysical method that allows one to estimate a 
value at early stages of the Universe evolution was orig- 
inally proposed by Savedoff (1956). Bahcall and Schmidt 
(1967) were the first to apply this method to fine-splitting 
doublets in quasar spectra. A recent update of this method 
has given the result: (—4.6 — > 4.2) x 10~ 14 yr -1 for up- 
per limit (95% C.L.) on a possible a variation (Cowie and 
Songaila 1995). A more stringent limit is presented in the 
paper by Varshalovich et al. (1996) and in this paper. 

We adhere the conventional belief that the quasar is 
an object of cosmological origin. Its continual spectrum 
was formed at an epoch corresponding to the redshift z 
of main emission details specified by relationship A Q b s = 

Alab(f + z). 
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Quasar spectra show the absorption resonance lines of 
the ions C IV, Mg II, Si IV, and others, corresponding to 
the S1/2 — ► -P3/2 (Ai) and Si/ 2 — ► P1/2 (A2) transitions 
(see Fig. 1). The difference between Ai and A2 is due to 
the fine splitting between the P 3 / 2 and P\/2 energy levels. 
The relative magnitude of the fine splitting of the corre- 
sponding resonance lines is approximately proportional to 
the square of a, 



SX 
T 



(3) 



where <5A = A 2 — Ai and A = (A 2 + Ai)/2. Thus, the ra- 
tio (5X/X) z /(SX/X)o is equal to (a z /a ) 2 , where the sub- 
scripts z and denote the fine-splitting doublet in a quasar 
spectrum and the laboratory value, respectively. Hence, 
the relative change in a can be approximately written as 



Aa \ 
a ) i 



{5X/X) Z 
(SX/X)o 



(4) 



provided that (Act/ a) is small. 

Thus, by measuring Ai and A 2 in an absorption sys- 
tem corresponding to the redshift z and comparing the 
derived and laboratory values, one can directly estimate 
the difference between a at the epoch z and the present 
value. 

It is to be noted that this method is the only one which 
provides direct values of Aa/a for different space-time 
points of the Universe. Thus, using this method we can 
study not only possible deviations of the fine-structure 
constant from its present value but also its values in space 
regions which were causally disconnected at earlier evolu- 
tionary stages of the Universe (e.g. Tubbs & Wolfe 1980; 
Varshalovich & Potekhin 1995). 

Tubbs & Wolfe (1980) (see also references therein) have 
used a coincidence of redshifts of optical resonant lines 
of ions with redshifts of the hydrogen 21 cm radio lines 
in distant absorption systems to derive an upper limit 
on the combination a 2 g p m c /m p at z ~ 0.4 — 1.8. Here, 
g p , m c and m p are the proton gyromagnetic factor and 
the masses of electron and proton, respectively. Recently, 
Drinkwater et al. (1998) used a similar method of com- 
parison of the redshifts of the HI (21 cm) and molecular 
radio lines at redshift z — 0.68 in two quasars to place a 
new strongest restriction on fractional variation of g p a 2 at 
the level <~ 10~ 15 yr _1 , from which they concluded that 
\a/a\ < 5 x 10~ 16 yr -1 . Two drawbacks are inherent to 
this method of comparison: (i) It does not give direct re- 
striction on variation of a, g p , or electron-to-nucleon mass 
ratios, but only on their combination, (ii) One cannot 
be sure that the hydrogen and molecular absorption lines 
originate from the same cloud along the respective lines 
of sight. Since the analysed line profiles were complex and 
required decomposition into several contours, the perfect 
juxtaposition of the hydrogen and molecular lines (hence 
the seemingly extra high accuracy) might result from acci- 
dental coincidence of redshifts of different components of 



the line profiles. Only analysis of the ratios of wavelengths 
of the same ion species (in particular, the doublet ratios 
mentioned above) is free of such ambiguity. 

Detailed discussion of possible sources of systemati- 
cal and statystical errors of the method based on Eq. (4) 
has been given elsewhere (Potekhin & Varshalovich 1994). 
The most significant source of possible systematical error 
turned out to be the uncertainty in the laboratory wave- 
lengths Ao (see below). Some errors which have a system- 
atic character for one selected absorption system (e.g., the 
shifts of estimated line centers resulting from occasional 
unidentified blending of partially saturated lines) become 
random (statistical) for a sample of unrelated absorption 
systems. For this reason, Potekhin & Varshalovich have 
argued that one should not confide in the errorbars esti- 
mated for individual absorption systems, but derive the er- 
ror from the actual scatter of the data. For example, Webb 
et al. (1998), using individual estimates of statistical er- 
rors of several fine-structure absorption systems, reported 
an unprecedented accuracy and claimed an evidence for 
nonmonotonic a variation: Aa/a = (— 2.64±0.35) x 10~ 5 
at 1.0 < z < 1.6. However, an independent statistical 
treatment of the data presented in their paper, disregard- 
ing the errorbars but using the actual scatter of the data 
to estimate the confidence level, yields (in the same range 
of redshift) Aa/a = (-3.0 ± 1.0) x 10~ 5 ± er syst , indi- 
cating that the individual statistical errors reported by 
the authors significantly underestimate the actual statis- 
tical uncertainty. Here, the systematic error cr sys t ~ 10 -4 
is mainly due to the uncertainty in the laboratory wave- 
lengths [see Eq. (5) below]. 

In this paper we apply the described method to a sam- 
ple of 20 absorption systems at 2.0 < z < 3.6, selected 
according to the criteria formulated in the next section. 

4. Criteria of data selection 

Of the available alkali-like doublets observed in quasar 
spectra such as C IV, Mg II, Al III, O VI, N V, Si IV, and 
others, we have selected for our analysis the Si IV line 
doublet, because it has the greatest ratio SX/X = 6.45 x 
10~ 3 , allowing this ratio to be measured most accurately. 
The abundance of silicon and its ionization state, as a rule, 
is such that the Si IV doublet lines occur on a linear part 
of the growth curve, which simplifies determination of the 
central wavelength of each line. 

The laboratory values of the Si IV doublet wavelengths 
(A01 = 1393.755 A and A 02 = 1402.769 A, according 
to Striganov & Odintsova 1982; A i = 1393.755 A and 
A 02 = 1402.770 A, according to Morton et al. 1988) are 
known with an uncertainty a\ i=a 1 mA. This uncertainty 
can introduce an appreciable systematic error in the de- 
termination of (Aa/a) z from Eq. (4). In the case of the 
Si IV doublet, this error can be estimated as 



a syst (Aa/a) w a x J(5XV2) « 8 x 10" 



(5) 
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This value of the systematic error for Si IV is the smallest 
among all alkali-like ion species listed above. 

Available observational data on Si IV doublets in 
quasar spectra have been selected taking into account the 
following criteria: 

1. High resolution FWHM < 25 km/s. 

2. Wavelength of each doublet component is measured in- 
dependently (without including the a priori informa- 
tion on SX/X). 

3. Any obvious blending of each component of the dou- 
blet is absent. 

4. Equivalent width (Wj.^) ratio of the doublet compo- 
nents is limited to 1 < (W x ± 2a Wl ) / (W 2 ± 2a W2 ) < 2 
(in conformity with the oscillator strength ratio). 

5. The lines are detected at the level higher than 5a, 
where a is a noise level. 

6. A dispersion curve for calibration is carefully deter- 
mined - e.g., using the technique described in the next 
section. 

Unfortunately, many of recent data found in literature 
do not often satisfy the necessary conditions 2, 3, and 5. 

5. Observations and data reduction 

Suitable observational data have been presented by Pe- 
titjean et al. (1994), Cowie & Songaila (1995), and Var- 
shalovich et al. (1996). To give an idea about basic ele- 
ments of observations and data reduction we describe here 
the main stage of the work by Varshalovich et al. (1996). 

Observations were carried out in 1994 using the Main 
Stellar Spectrograph of the 6-m Telescope. A Schmidt 
camera (F:2.3), reconstructed to operate with a CCD ar- 
ray, was employed. In combination with a 600 lines mm -1 
diffraction grating operating in the second order, this 
camera yielded a linear reciprocal dispersion of 14 A 
mm" - 1 (0.24 A pixel" 1 ). Since this study is metrological 
in essence, special attention was paid to wavelength cali- 
bration and check measurements. A Th+Ar lamp and the 
Atlas of the Th+Ar spectrum constructed from echelle 
data (D'Odorico et al. 1987), which gave vacuum wave- 
lengths A to within 0.001 A, were employed for calibra- 
tion. Ten to twenty reference lines were used to construct 
the dispersion curve in the form of a Chebyshev poly- 
nomial (A = aoTo + a\T\ + C12T2 + where T n is the 
Chebyshev polynomial of nth order) for each of the spec- 
tral segments in question. The nonlinear terms (a 2 ) of the 
dispersion curve were determined especially carefully, be- 
cause the accuracy of the measured ratio SX/X was partic- 
ularly sensitive to the nonlinearity of the scale, contrary 
to the absolute calibration (i.e. ao). For the same reason, 
the spectrograph was adjusted for each of the Si IV dou- 
blets in such way that both doublet components, with a 
separation of about 36 A for z ~ 3, were in the middle 
of the spectral segment corresponding to one diffraction 
order, where the nonlinear distortions of the scale were at 
a minimum, while the sensitivity was at a maximum. 
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Fig. 2. The spectrum of HS 1946+76: comparison of the pro- 
files of the Si IV doublet absorption lines with z a b s = 3.049 
[the Ao = 1393 A (a) and 1402 A (b) lines in our study and 
the Ao = 1402 A line (c) in Fan and Tytler (1994)] 

The wavelengths Ai and A2 were determined by gau- 
sian decomposition of the corresponding absorption lines 
in the observed quasar spectra. Portions of the spectrum 
of HS 1946+76 (z cm = 3.05), one of the brightest QSOs in 
the sky, measured by Varshalovich et al. (1996), are shown 
in Figs. 1 and 2. Despite the absorption doublet shown in 
Fig. 1 consists of two components, it still satisfies the cri- 
teria formulated in Sect. 4, including item 3. In fact, we 
have checked that even total neglect of the second com- 
ponent of each line alters the estimate of Aa/a for this 
absorption system by less than 2 x 10~ 5 , which is well be- 
low the overall statistical uncertainty reported in Eq. (6) 
below. On the other hand, the portion of spectrum shown 
in Fig. 2 requires caution, since the spectral components 
are close and partially saturated. 

In their recent analysis of the spectrum of HS 1946+76, 
Fan and Tytler (1994) have detected six absorption sys- 
tems, two of which contain Si IV lines. However, they pre- 
sented only one of the doublet components in each of the 
doublets. Therefore, their high-resolution results cannot 
be used for measuring Aa. In Fig. 2, the profiles of our 
measured spectral lines are compared with the profile ob- 
tained by Fan & Tytler with resolution of ~ 10 km s . 
The comparison confirms that the resolution achieved in 
our work is equally sufficient for the profile analysis. 

6. Results 

6.1. Bounds on a 

Results of the analysis of the Si IV fine-splitting doublet 
lines are presented in Table 1. The third column shows a 
deviation of the a value calculated according to Eq. (4) 
for a single doublet. The value (<5A/A) = 0.0064473 in 
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Eq. (4) has been adopted from Morton ct al. (1988). From 
the data listed in Table 1, (Aa/a) is estimated by stan- 
dard least squares. According to the arguments given 
in Sect. 3, we did not involve estimates of the individ- 
ual wavelength uncertainties. For example, for the profile 
shown in Fig. 1, the uncertainty of the gaussian decom- 
position yields a\ ~ 0.007 A, and the estimate a\ as a 
function of resolution and signal-to-noise ratio, following 
Young et al. (1979), yields even smaller a\ <~ 0.005 A, 
while the actual uncertainty may be larger. Estimating the 
statistical error from the scatter of Aa/a measurements 
and taking into account the systematic error in Eq. (5), 
we arrive at the final estimate (for z interval 2.0 — 3.5): 

(Aa/a) = (-3.3 ± 6.5 [stat] ± 8.0 [syst]) x 10~ 5 . (6) 

The corresponding bound at the 95% significance level is 



\Aa/a\ < e = 2.3 x 10" 



(7) 



Since we have considered a sample of 20 absorption 
systems, the statistical (observational) error in Eq. (6) is 
~ \/20 ~ 4.5 times smaller than a typical error of a sin- 
gle measurement. On the other hand, the statistical error 
and the systematic error due to the above-mentioned un- 
certainty of the laboratory wavelengths are of the same 
order of magnitude, hence any significant further improve- 
ment of the upper limit (7) would require an improvement 
of not only observational techniques but also the accuracy 
of the laboratory measurements. 

According to the standard cosmological model with 
the cosmological constant A = 0, the time elapsed since 
the formation of the absorption spectrum with redshift z 
is 



Mi 
Mi 



(1 + z)- 3 / 2 ] 

(i+^n 



(fi = l), 
(fi« 1), 



(8) 
(9) 



where f2 is the mean-to-critical density ratio. With the 
present age of the Universe equal to (14 ± 3) x 10 9 yr, we 
see that the redshifts z = 2.8 ± 0.7 presented in Table 1 
correspond to the elapsed time ~ 12 ± 4 Gyr. 

Thus, the upper limit on the time-averaged rate of 
change of a (with most likely values z = 2.9 and t = 
14 x 10 9 yr) at the 2a level is 



a/a 



< 1.9 x 10~ 14 yr 



(10) 



And the most conservative upper limit (z 
11 x 10 9 yr) at the 2a level is 



2.1 and t 



a/c 



< 3.1 x 10 



-14 



yr 



(11) 



Table 1. Variation of a values estimated from redshifted 
Si IV doublet line splitting according to Eq. (4) 



Quasar 


2 


/A t \ -1 f-^ A. 

(Aa/a) z ,W 


Ref. 


HS 1946+76 


3.050079 


1.58 


1 


HS 1946+76 


3.049312 


0.34 


1 


HS 1946+76 


2.843357 


0.59 


1 


S4 0636+68 


2.904528 


1.37 


1 


S5 0014+81 


2.801356 


-1.80 


1 


S5 0014+81 


2.800840 


-1.70 


1 


S5 0014+81 


2.800030 


1.11 


1 


PKS 0424-13 


2.100027 


-4.51 


2 


Q 0450-13 


2.230199 


-1.48 


2 


Q 0450-13 


2.104986 


0.02 


2 


Q 0450-13 


2.066646 


1.03 


2 


Q 0302-00 


2.785 


2.07 


3 


PKS 0528-25 


2.813 


1.29 


3 


PKS 0528-25 


2.810 


1.03 


3 


PKS 0528-25 


2.672 


-5.43 


3 


Q 1206+12 


3.021 


-1.29 


3 


PKS 2000-33 


3.551 


-3.88 


3 


PKS 2000-33 


3.548 


2.85 


3 


PKS 2000-33 


3.332 


5.95 


3 


PKS 2000-33 


3.191 


-5.69 


3 



References: [1] Varshalovich et al. (1996); [2] Petitjean et al. 
(1994); [3] Cowie and Songaila (1995). 



6.2. Consequences for theoretical models 

The results obtained above make it possible to select a 
number of theoretical models for dependence of a on cos- 
mological time t (the time elapsed since the Universe cre- 
ation). Consider some of them: 
a. The hypothesis of the logarithmic dependence 



In 



t 



87TT Planck 



(12) 



Here, Tpi anc k = \/h~G/c 5 = 5.4 x 10~ 44 sec is the Planck 
time. Teller (1948) suggested this dependence on the ba- 
sis of the striking concurrence of the value a -1 = 137.036 
and value ln(io/87rTpi anc k) ~ 137 (note that the form of 
the dependence and t value is presented in the modern in- 
terpretation). Later Dyson (1972) showed that it could be 
supported by considerations that follow from the method 
of renormalization in the QED. In the standard cosmolog- 
ical model with A = and fi = 1, this hypothesis gives 
rise to the following dependence 

a a 



a z = 



l-|a ln(l + z) 



(13) 



For z — 2.8, this yields a z = 1.015 ao, in obvious contra- 
diction with inequality (7). Thus, the logarithmic depen- 
dence can be excluded. Note that a more general depen- 
dence a^ 1 = C\n(t/T) with arbitrary C and r > Tpi anc k 
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has been ruled out already by Varshalovich & Potekhin 
(1995). 

b. The power-law dependence (for standard model) 



(14) 



in the standard flat Universe this leads to the dependence 
a = a (l + z)~ Zn / 2 , from which, using Eq. (7), we obtain 
the bound 



\n\ < 2e/[31n(l + z)} = 1.1 x 10" 



(15) 



c. The Kaluza-Klein-like models often assume that com- 
pactification was occurring at the very early stages of 
Universe evolution and now the additional D dimensions 
of the subspace have radius of curvature R ~ Zpi an ck ~~ 



\/hG/i 



1.6 x 10 33 cm. Thus, determining of varia- 



tions at such level is impossible at present. On the other 
hand, there are theories in which the process of the com- 
pactification can develop even now. The information about 
R for such theories may be obtained from measurements 
of the a variation. According to Eq. (2), the restriction on 
the changing of the subspace scaling factor is 



R 2 a 



(16) 



d. The recent version of the dilaton evolution proposed 
by Damour & Polyakov (1994) in the frame of the string 
model has provided an expression connecting the variation 
of the fundamental constants: 



G 
G 



Ok 



a 
a 



3332. 5fc- 



a 



(17) 



Here, k is a parameter of the theory, which fixes a type 
and speed of the evolution of the scalar field (i.e. dilaton 
field), 6 = 2(ln(A s /mc 2 )) 2 is a numerical coefficient, A s ~ 
5 x 10 17 GeV is a string mass scale, and m = 1.661 x 
10~ 24 g is the atomic mass unit (mc 2 = 931.5 MeV). The 
parameter k regulates the character of the time variation: 
G and a oscillate for k > 1, whereas for k <C 1, these 
constants change monotonously. For the models with k > 
1 discussed by Damour & Polyakov (1994), a measured 
bound on G/G enables one to obtain a bound on a/a. 
For models with k <C 1, we can obtain the bound on G/G 

from Eq. (11): G/G < 3332.5 \a/a\ < 1 x lO" 10 yr" 1 . 



7. Conclusions 

We have formulated criteria for selection of high-quality 
doublet lines in quasar spectra, most suitable for the anal- 
ysis of the large-scale space-time variation of the fine- 
structure constant. We have also described the calibration 
technique capable to provide the necessary metrological 
quality. 

Using the formulated criteria and described technique, 
we have performed the Si IV doublet analysis of quasar 



spectra, that has enabled us to set an upper bound on 
the rate of the possible cosmological time variation of the 
fine-structure constant. This limit is much more model- 
independent and hence more robust than those provided 
by the local tests, including the analysis of the Oklo phe- 
nomenon. 

One may anticipate that with development of observa- 
tional spectroscopy the portion of material which satisfies 
the formulated criteria will rapidly increase in the near 
future. However, as shows the discussion in Sect. 6.1, this 
will not allow one to further tighten the bound on a/ a 
until an improvement of the laboratory wavelength accu- 
racy. On the other hand, in frames of the hypothesis that 
a does not change at the present stage of the cosmological 
evolution, the progress of the observational spectroscopy 
offers an exciting opportunity to improve the accuracy of 
determination of the fine splitting 6X/X against that avail- 
able in the laboratory experiments. 
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